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ABSTRACT 


A system which is either working or under repair is modeled as an alternating re- 
newal process. Data is in the form of a finite number of independent lifetimes and repair 
times. Several semi-parametric estimators of the probability that the system is up at a 
finite time t are studied via simulation. The estimators are based on an exponential ap- 
proximation to the true system availability at time t and use empirical Laplace trans- 


forms of the lifetimes and repair times. 
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I. INTRODUCTION 


Assume a new system has been installed at t=0 and 1s up. The system is subject to 
failure and repair. An initial simple model for this situation is an alternating renewal 
process. 

Let U, be the /* time between failures and D, be the i* repair time. Assume {U,} are 
independent identically distributed random variables with distribution F and {D,} are 
independent identically distributed random variables independent of {U,} with distrib- 
ion G. Let 


X,= | if the system is up at time t 
0 otherwise. 
A renewal argument yields 


A(t) = P(X, =11X% =1} 


=1- FU + [Pod )40-y) (1.1) 


Where F * G denotes the convolution of F and G. It can be shown that the long run 


proportion of time the system is up is 


A(co) = lim A(t) = ae 
as E(U,) + E(D,) 
If F 1s exponential with mean /-! and G is exponential with mean pu“! the equation 
(1.1) has the solution 


ale | 
NQ) re eas 
rues ot 


However, for most distributions F and G, equation (1.1) has a very complicated closed 
form solution. 


Computing the Laplace transform of equation (1.1) 


A(s) = | aa | e-"'F(t)dt + \ "(PGMA AU~ 9) 


0 


+ fl- F(s)] + F(s)G(s)A(s) 


aoe | e~"G(dt) . 
0 
Solving the equation results in 


l A 
A apele etS)3 
A> eee 
1 — F(s)G(s) 
The following approximation to A(t) 1s proposed 
A,(t) = A(oo) + (1 — A(oo))e*" 


for some f . One choice for f is to set 


|, se"'A(tdt = |. se A(tjdt 
and solve for f ats= +. More specifically let 


A(t) — Aloo 
Oe ae 


and set 


(1.2) 


(1.5) 


and 
ae l mA, (1.6) 
P(s) 
where 
Gj = (1.7) 
1 — F(s)G(s) 
ape 5 = = 


The choice of £ can be “tuned” to t by use of a different weight function. For ex- 


ample, f can be chosen so that 


| _ w()A Anat = | ~ w(NA(at (1.8) 


where w(t) is the density function of the median of 5 independent exponential random 
variables each having mean + and s 1s chosen so that the expected value of the median 
equals t. The resulting expression for f involves the Laplace transforms of F and G. 
Other approximations can be found in Gaver and Jacobs[Ref. 1]. 

Suppose that data are available on the successive up times and repair times of the 
SVStem; Sav uj, UW, ...,u, and d,, d@,, ....d,. A problem of practical interest 1s to use these data 
to estimate the probability the system will be up at a finite time t, often referred to as 
the availability of the system at time t. One estimate that is often used in practice is the 
ratio 

ud 


(u +d) 


where 


my 
ll 
x|— 


ui; 


and 


a ee | 

d=— Z an 
This is an estimator of the long-run proportion of time the system is up. This estimator 
can be quite biased for small times, t. 

In this thesis estimators of the availability of the system at time t are investigated 
bv simulation. Three of the estimators are based on the exponential approximation 
(equation (1.3)) and use empirical Laplace transforms to estimate the Lapiace trans- 
forms of F and G. Another estimator that is investigated is one which simulates the al- 
ternating renewal process using the observed data and a bootstrap-like sampling scheme. 

The specific estimators are described in Chapter 2. Chapter 3 contains details of the 


simulation experiment and results. Conclusions from the study are given in Chapter 4. 


II. NATURE OF THE PROBLEM 


A. PROBLEM 

A system subject to failure and repair is modeled by an alternating renewal process. 
The problem is to estimate the system availability at time t from a sample of lifetimes, 
U,, U;,...,u, and repair times d,,d,, ...,d,. In this chapter four estimators will be described. 
The first is a computationally intensive bootstrap estimator. The last three estimators 


are based on the exponential approximation (equation(1.3)). 


1. Empirical Estimator 

This estimator simulates the alternating renewal process using the observed data 
U,, Uz, ...,u, and d,, d,, ...,d@,. Specifically, a lifetime is drawn at random with replacement 
from {u,}. If the lifetime is greater than t, the system is said to be up. If the lifetime 1s less 
than t, a repair time is drawn at random with replacement from {d} and added to the 
lifetime. If the resulting sum is greater than t, the system is said to be down. If the sum 
is less than t another lifetime 1s drawn at random with replacement from {u,} and added 
to the sum. If the resulting sum is greater than t, the system 1s said to be up. If it is less 
than t, a repair time is drawn at random with replacement from {d}, etc. This procedure 
is repeated ., = 100 times and the fraction of times the system is up 1s computed. This 
is the empirical estimate of availability Asnp(t). Confidence intervals are obtained for this 
estimate by repeating the whole procedure AN’, times. The order statistics of the resulting 
\, availability estimates are used to obtain confidence intervals. Specifically the 
ey ).V,% and (1 — a ).V," order statistic are used to compute two-sided (1 - a)%o confi- 
dence interval, and the aN,” order statistic and | are used to compute one-sided (1 - a 
)% confidence interval. The procedure described is a Bootstrap, as discussed extensively 
by Efron [Refs. 2,3]. 


2. Exponential estimator 
The exponential estimator is based on the exponential approximation 
(equation(1.3)) with weighting function the exponential density resulting in equation 
(1.6). The estimator is obtained by replacing A(oo), F(s) and G(s) by the following esti- 


mates. 


i 
Tae 


Pom, 





Hine 


3 


F(s) is estimated by the empirical Laplace transform 


ou(s) = a ye 
i=] 


G(s) is estimated by 


and s 1s taken to be rs. The availability at time t 1s estimated by 


NO Gj es ee 


where 
eal 
A | Lm 
b= a I 
oe), 
with 
| + 4(+)— Aloo) 
A l 
ora are - 
~- (1 — A(oo)) 
and 
1-¢y+) 
Je t 
t rea 


1-¢4+)6(+) | 


A Taylor expansion of equation (1.6) for small s vields 


(2.1) 


(2.2) 


(2.3) 


(2.4) 


(2.5) 


(2.6) 


3) 


PING ELC |e EV )} 
a a as s—>Q (2.8) 
PC ee — Lene sy) 
where 


If B.> 100 then it is set equal to the sample version of equation (2.8); that 1s 
ELU,] is estimated by u, E[C,] is estimated by 7 + d, Var[C,] is estimated by 


ae | _ =) ] oR? 
oO. = N ae l Y u) as N aa l Ya da) (2.9) 
i= i= 


and E[C}] is estimated by 62+ (u4+d/). 
If £. iS negative it 1s set equal to = + 





g|— 


3. The Simplified Exponential Estimator 


The simplified exponential estimator has the form 


Win) = Alcs) + (1 — Aleo)e7 8 (2.10) 


> 


's larger than 100 


§ 


where B. is as computed in equation (2.5). It differs from B. in that, 
as before. 


x1] 


are not recomputed. However, if B. <0, than it is set equal to + + 
4. Cubic Estimator (Median-of-5 Exponential Approximation) 


Let Z have the distribution of the median of 5 exponentials each having mean 


l — 
|. Its density function is 


w= PZ <0) (2.11) 
4 
oe s(3)e = or lsrewi ea.” 
a 30r[e >" — pemee ate egy 


The expected value of Z is 


ee eee ik 
AEA ea ap a) en) 


set 


pclae 
lar (2.13) 


to ensure E[Z] = t. The cubic estimator is of the form 
GV AGS 4h = Zeal (2.14) 


where B. is the solution of the sample version of equation (1.8) with weight function w(t) 
the density function of the median-of-5 exponentials (equation (2.11) ). More specif- 


ically, 
if 
_ AW) = Aloo) _ -pr 
P(j= T= eS ~e (Dai) 
then setting 
[. #9Pldt = | w(te dt (2.16) 
0 


results in a cubic equation for # since 


[. mPa = 30rl | e"P(dt — 2 | e" P(ndt + | e' P(ndr] (2.17) 
0 0 0 


= 30r[ P(3r) — 2P(4r) + P(Sr)] 


and 


0 


| w(te edt = 30rf | ee = | Cn ae agi | ea (2.18) 
0 0 0 


2r° 


= 301. FF pyar + Byer + B) 


ie 


Setting equations (2.17) and (2.18) equal results after some algebraic manipulation in the 


cubic equation for f 


2 
fb) = B? + 12rB? + 4777p + 607? “> __—_ = 0. (2.19) 
P(3r) — 2P(4r) + P(5r) 
The first derivative of f( B ) 1s 
Sf? (B) = 38? + 24rB + 47r’ . (2.20) 


To solve the cubic equation (2.19) the quadratic equation f’ (f) =0 is solved for B to 


find the B ’s associated with the minimum and maximum. The resulting f ’s are 
AL 
B,=(-4-“yv (2.21) 


AV 
(44> —)r 


Note that 8, and £, are alway negative. Hence f({ B ) = O has at most one positive sol- 


ution or all negative solutions. If f(0) > 0, there is no positive solution and 


sj|— 


, aa (2.22) 


si] 


If f(0) < 0, there is a positive solution. To find f, for f{ B, ) = 0, the Newton search 
method[Ref. 4] is used. The availability at time t is estimated by 


Aco) i Atea)en | (2.23) 
B. CONFIDENCE INTERVALS 


Suppose X is a random variable whose probability law depends on an unknown 


parameter @. Given a random sample of X ; x, x, %;,...,x,, tWO Statistics lower 


(L = L(x) = L(%, %,...,x,)) and upper (U = U(x) = U(x, x,,...,x,)) form a 100(1- a )% 
confidence interval for 6, ifin repeated random samples L<@<U_ 100(1 — a) % of the 
time. 

To obtain confidence intervals for the exponential approximation point estimates for 
P{system up at time t], the Bootstrap Estimation method is used. Efron(1979) intro- 
duced the Bootstrap method to estimate the distribution of a statistic computed from 
observations [Refs. 2,3]. In this thesis the Bootstrap is implemented as follow: 

Given a sample of n lifetimes u,, %,...,u, and n repair times d,, d,,...,d,, sample at 
random with replacement n times from {u,} and n times from {d,}. This constitutes one 
bootstrap sample. Compute the estimate of interest 4,(t), A,(t) , or A,(t) from the boot- 
strap sample. Generate another independent bootstrap sample etc. until there are N, 
estimates. Order the N, estimates. The (=)NP and (1 ee order statistics of the 
estimates form the two-sided (1 - «)% confidence interval of A(t). The one-sided (1 - « 
)% confidence interval has a lower interval point of the «N," order statistic and 1 as the 


upper interval point. The number of bootstrap replications is NV, = 100. 


Ii. ANALYSIS OF THE PROBLEM 


A. SIMULATION 

A Fortran program is written to generate and analyze the data for this problem. All 
simulations are carried out on an IBM3033AP computer at the Naval Postgraduate 
School using the LLRANDOM II random number generating package [Ref. 5]. 

The system lifetimes are independent identically distributed exponential with mean 


A-' = ], The system repair times have one of three distributions in this study: 
1. Exponential with mean > 
P{R>nj=e" psn): 
2. Gamma 
P{R>thh=e7' + 41e" t > 0; short-tailed; increasing hazard; 
3. Mixed Exponential 
P{R>t= 0.9e7 7 + 0.le72! t>0; long-tailed; decreasing hazard. 


Table 25(Appendix B) contains the true availability A(t) at various times t for each of 
these models. Tables 26 through 31(Appendix C) contain the mean, and square root of 
mean square error(SRMSE) of the estimated availability Am), A(t), A(t) and A(t) 


at Various times t for each of these models; i.e. the mean is 





where A,(t), 1s the point estimate at time t in the k* super-replication. 
The simulation has N, = 500 super-replications. For each super-replication N = 50 
or 100 lifetimes and repair times are generated. Using this sample the Empirical Esti- 


mate, Exponential Est:mate, Simplified Exponential Estimate, and Cubic Estimate are 
computed. A, = 100 bootstrap replications are drawn and the estimates are recomputed. 
The 90 % two sided confidence intervals and 90 % one sided confidence intervals are 
computed using NV, = 100 bootstrap estimates. Recorded in Tables 1, 3, 5, 7, 9, 11, 13, 
15, 17, 19, 21 and 23 are the number of replications for which the confidence intervals 
cover the respective true probabilities; also recorded are the number of confidence in- 
tervals which are too low (true P[system up at time t] > upper bound) and too high (true 
P[system up at time t] < lower bound). The average length of the confidence interval is 
computed as well as the standard deviation of the lengths. The standard deviation 1s 
computed by subtracting the mean length from each length, squaring the results, sum- 
ming over the 500 lengths and dividing by 499, and finally taking the square root of the 


result. 


B. ANALYSIS & RESULTS 

In this section results will be reported for the simulation experiment. Some true 
values of P[system up at time t] for these models can be found in Appendix A. The 
number of bootstrap replication is set at 100. The simulation 1s replicated 500 times. For 
each replication, 90 % two sided confidence intervals and 90 % one sided confidence 
intervals of P[system up at time t] are computed using each procedure of Chapter 2. The 
times considered are t = 0.2, 0.5, 1.0, 2.0 and 3.0. The sample size N is set at 50 and 
100. 

The confidence interval procedures for P[svstem up at time t] use the Bootstrap for 
all estimators. For each procedure, the number of intervals covering the true value of 
P[svstem up at time t] 1s recorded as well as the number of intervals that are too high 
or too low. These results are reported in tables J], 3, 5, 7, 9 and 11 for two sided confi- 
dence intervals, and tables 13, 15, 17, 19, 21 and 23 for one sided confidence intervals. 
Next to each coverage count 1s given the corresponding coverage proportion in paren- 
thesis. If a (1- « ) % confidence interval procedure is performing well, then this interval 
should cover about (l- a ) % of the time. In tables 2, 4, 6, 8, 10 and 12, the average 
length of the two sided confidence intervals for P[svstem up at time t] are recorded, and 
in tables 14, 16, 18, 20, 22 and 24, the average length of the one sided confidence inter- 
vals for P[svystem up at time t] are recorded. The estimated standard deviation of the 
length is below the average length. If an estimator performs well, its confidence interval 


should not only have the correct coverage rate but also a small average length. The 


simulation results recorded in tables 1 through 4 and 13 through 16 are for the Expo- 
nential repair time. Tables 5 through 8 and 17 through 20 for the Gamma repair time. 
Tables 9 through 12 and 21 through 24 for the Mixed exponential repair time. Coverage 
results for sample size N = 50 are presented in table 1, 5, 9 13, 17 and 21. Sample mean 
and standard deviation of the confidence interval lengths for sample size N = SO are 
presented in table 2, 6, 10, 14, 18 and 22. Coverage results for sample size N = 100 are 
presented in table 3, 7,1, 15, 19 and 23. Sample mean and standard deviation of the 
confidence interval lengths for sample size N = 100 are presented in table 4, 8, 12, 16, 
20 and 24. 

To assess the coverage results, recall that for a binomial random variable with pa- 
rameters 500 and p = 0.9, a 95 % confidence interval for p is (0.874, 0.926] or for the 
number of successes 1s [436.9, 463]. 

The Empirical Estimator tends to undercover for the Exponential and Gamma re- 
pair time cases for sample size N = 50 and overcover for all repair time distributions for 
sample size N = 100. For all repair time distributions, the coverage rate and the con- 
fidence interval width of the Empirical Estimator tend to increase according to increasing 
time t. The Empirical Estimator tends to have the largest confidence interval for all re- 
pair time distributions. The simulation results for time t = 2.0, 3.0 at sample size N = 
100 are not available because of the very large computational requirements of the pro- 
cedure. 

The two-sided 90 % confidence intervals for the Exponential Estimator tend to have 
the correct coverage for all repair time distributions at sample size N = 50: that is, the 
number that cover are within the Binomial 90 % confidence interval [437, 463]. In- 
creasing the sample size to N = 100 results in under-coverage, suggesting the estimator 
is biased. For all repair time distributions, two-sided confidence intervals tend to have 
better coverage than one-sided confidence intervals at sample size N = SO, and the 
one-sided confidence intervals tend to have better coverage than the two-sided confi- 
dence intervals at sample size N = 100. 

The Simplified Exponential Estimator tends to undercover for all repair time dis- 
tributions. The two-sided coverage rate of the Simplified Exponential Estimator for 
Mixed exponential repair time tends to decrease and the confidence interval width tends 
to increase according to increasing time t at sample size N = 50. The Simplified Expo- 


nential Estimator tends to have the smallest confidence intervals for Exponential and 
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Gamma repair time. The best coverage for these confidence intervals is for the Expo- 
nential repair time at sample size N = 100. 

The two-sided confidence intervals of the Cubic Estimator tend to undercover for 
all repair times. The coverage rate of the Cubic Estimator for exponential repair time 
tends to increase and for Mixed exponential repair time tends to decrease according to 
increasing time t at sample size N = 100. The confidence interval width tends to in- 
crease according to increasing time t. 

For two-sided confidence intervals, the Exponential Estimator has the best coverage 
of the three estimators based on the exponential approximation. The Cubic Estimator 
confidence intervals tend to have slightly better coverage than the Simplified Exponen- 
tial intervals. For all estimators, the average length of the confidence intervals is the 
largest for the Mixed exponential repair time. 

The one-sided confidence interval coverage is about the same for all three estimators 
based on the exponential approximation. The results of the tables in Appendix B indi- 
cate that the means of all estimators tend to be close to the true availabilities, A(t), for 
all repair time distributions. The root-mean-square error of the Exponential Estimator 
is the smallest for most cases. Increasing the sample size from N = 50 to N = 100 re- 
sults in more accurate estimates and decreases the root-mean-square error, as would be 
anticipated. 


IV. CONCLUSIONS 


This thesis considers the problem of estimating the availability of a system at a finite 
time t. Simulation is used to study four estimation procedures and associated confidence 
interval procedures. The basic model is the alternating renewal process; the successive 
lifetimes are independent identically distributed and the successive repair times are in- 
dependent identically distributed. Three repair time distribution are used in the simu- 
lations: the exponential, the gamma, and the mixed exponential; the lifetime distribution 
1s exponential. 

Three of the estimators are based on an exponential approximation to the avail- 
ability at time t. The fourth estimator, called the Empirical Estimator, simulates the 
underlying alternating renewal process using a bootstrap-like sampling scheme. 

The following conclusions are drawn from the simulation experiment. 


1. The two-sided confidence intervals of the Exponential Estimator have the best 
coverage among the two-sided confidence intervals. 


2. When sample size is increased from N = 50 to N = 100, the two-sided confidence 
intervals tend to undercover. When a confidence interval does not cover, it tends 
to be too high indicating that the estimators based on the exponential approxi- 
mation may be biased towards overestimating the availability. 


3. For a sample size of N = 50, the two-sided confidence intervals for the Empirical 
Estimator tend to undercover. For N = 100, they tend to overcover. The one-sided 
confidence intervai coverage for the Empirical Estimator is better than its two-sided 
confidence interval coverage. The coverage of the confidence intervals of the Em- 
pirical Estimator may be improved by increasing the number of bootstrap-like 
replications, V, , used to compute the estimate. 


4. The computational effort in computing confidence intervals for the Empirical Es- 
timator is the greatest and may be prohibitive for larger times. 


5. The average two-sided confidence interval lengths are the largest for Gamma repair 
times and the smallest for Mixed exponential repair times. This behavior 1s due to 
the difference in the true system availability under the different assumptions. 


6. The Exponenual Estimator has the smallest root-mean-square error for the Expo- 
nential repair ttme and Gamma repair time cases of all the estimators. 


7. Of the four estimators, the Exponential Estimator appears to be the most promis- 
ing. 
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APPENDIX A. COVERAGE RATIO & LENGTH OF C.I TABLES 


Table 1. COVERAGE RATIO (EXPONENTIAL REPAIR TIME : TWO-SIDED, 
N = 50) 


. Eh Ges ree 
















3.0 






439(0.88) 442(0.88) 418(0.84) 420(0.84) 
24(0.05) 20(0.04) 47(0.09) 45(0.09) 


Table 2. LENGTH OF C.J (EXPONENTIAL REPAIR TIME : TWO-SIDED, 
N = 50) 


0.2 
0.1380 0.1354 O62 On! 

ee 
0.1476 0.1349 0.1324 0.134 

1.0 
0.1476 0.1363 OnsSd 0.1342 

2.0 
0.0153 0.0196 0.0186 0.0185 

3.0 
0.0143 0.0201 0.0185 0.0184 













Table 3. COVERAGE RATIO (EXPONENTIAL REPAIR TIME : TWO-SIDED, 
N = 100) 


2.0 
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Table 4. LENGTH OF C.I (EXPONENTIAL REPAIR TIME : TWO-SIDED, 
N = 100) 


0.1116 0.0866 0.0846 0.0877 
0.2 

0.0156 0.013 0.0127 0.012 
0.5 

0.0139 0.0105 0.010 0.0098 

0.1475 0.0974 0.0948 0.0962 
2.0 
3.0 













Table 5. COVERAGE RATIO (GAMMA REPAIR TIME : TWO-SIDED, 
N = 50) 


te 












3.0 






461(0.92) 450(0.90) 431(0.86) 431(0.86) 
17(0.03) 21(0.04) 43(0.09) 43(0.09) 


Table 6. LENGTH OF C.I (GAMMA REPAIR TIME : TWO-SIDED, N = 50) 


0.2 

0.1425 0.1328 0.1284 0.1352 
0.5 

0.0152 0.0159 0.0147 0.0167 
2.0 

0.149 0.1196 0.117] 0.1173 
3.0 
0.0144 0.0171 0.0159 0.0159 
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Table 7. COVERAGE RATIO (GAMMA REPAIR TIME : TWO-SIDED, 
N = 100) 


38(0.07)1 35(0.07) 35(0.07) 36(0.07) 
0.2 
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Table 8. LENGTH OF C.I (GAMMA REPAIR TIME : TWO-SIDED, N = 100) 


0.2 

0.1436 0.096 0.0937 0.0969 
» 0.5 

0.1485 0.0879 0.0857 0.0878 
2.0 
3.0 
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Table 9. COVERAGE RATIO (MIXED EXPONENTIAL REPAIR TIME : 
TWO-SIDED, N = 50) 


too low 25(0.05) 9(0.02) 
Z0 









too high 
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Table 10. LENGTH OF C.I (MIXED EXPONENTIAL REPAIR TIME 
TWO-SIDED, N = 50) 


0.0895 0.0727 0.0708 0.0715 

0.2 
0.0181 0.0153 0.0148 0.0147 
0.0997 0.0751 0.0732 0.0728 

0.5 | 

0.0157 0.016 0.0153 0.0147 
0.0149 0.0203 0.0197 0.0187 

2.0 
0.0152 0.0246 0.0244 0.0238 

3.0 
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Table 11. COVERAGE RATIO (MIXED EXPONENTIAL REPAIR TIME : 
TWO-SIDED, N = 100) 


29(0.06) 33(0.07) 34(0.07) 34(0.07) 
0.2 









eae [sew on ol 
feenm | | wean | om [ae 
3.0 






} 100 Iow 


Scorer | 426(0.85) 436(0.87) 436(0.87) 
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Table 12. LENGTH OF C.I (MIXED EXPONENTIAL REPAIR TIME 
TWO-SIDED, N = 100) 


0.0904 0.0541 0.0526 0.0539 
Or2 

0.0142 0.0078 0.0077 0.0078 
0.5 

0.1051 (0595 0.058 0.0573 

OOles 0.0125 0.0115 0.0109 
2.0 
BAU 













Table 13. COVERAGE RATIO (EXPONENTIAL REPAIR TIME: ONE-SIDED, 
= 50) 
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Table 14. LENGTH OF C.I (EXPONENTIAL REPAIR TIME : ONE-SIDED, 
N = 30) 


0.2 
0.0241 0.0281 0.0281 0.0267 
0.1585 0.1521 0.1511 0.1542 
0.5 
0.0219 0.0233 0.0224 0.0226 
240 
3.0 
0.0213 0.0253 0.0248 0.0247 













ZS 


Table 15. COVERAGE RATIO (EXPONENTIAL REPAIR TIME: ONE-SIDED, 
N = 100) 
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Table 16. LENGTH OF C.I (EXPONENTIAL REPAIR TIME : ONE-SIDED, 
N = 100) 


OsZay 0.0959 0.0959 0.1003 

0.2 
0.1585 Oven 0.1105 0.1126 

0.5 
0.1703 0.1104 0.1105 0.1123 

2.0 

3.0 
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Table 17. COVERAGE RATIO (GAMMA REPAIR TIME : ONE-SIDED, 
N = 50) 


estas a pec 
0.5 


cover 401(0.80) 405(0.81) 394(0.79) 406(0.81) 











too high 69(0.14) 81(0.16) 96(0.19) 95(0.19) 
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Table 18. LENGTH OF C.I (GAMMA REPAIR TIME : ONE-SIDED, N = 50) 


0.1302 0.1412 0.1429 0.1518 
0.2 

0.0234 0.0296 0.0289 0.0286 

0.1642 0.1518 0.151 0.1575 
0.5 

0.0215 0.0212 0.022 0.0234 

0.0226 0.0201 0.0199 0.0203 

0.1731 0.135 0.1335 0.1341 
2.0 

0.1729 0.1344 0.1331 0.1333 
3.0 
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Table 19. COVERAGE RATIO (GAMMA REPAIR TIME : ONE-SIDED, 
N = 100) 
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Table 20. LENGTH OF C.I (GAMMA REPAIR TIME : ONE-SIDED, N = 100) 


0.2 
0.0206 0.0165 0.0165 0.0166 

0.5 
0.0204 0.0139 0.0137 0.0149 
0.0197 0.0126 0.0128 0.0133 

2.0 

3.0 


















_ 


Table 21. COVERAGE RATIO (MIXED EXPONENTIAL REPAIR TIME : 
ONE-SIDED, N = 50) 
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Table 22. LENGTH OF C.I (MIXED EXPONENTIAL REPAIR TIME 
ONE-SIDED, N = 50) 


0.2 
0.0222 0.0166 0.0166 0.0161 

0.5 
0.0194 0.0217 0.0215 0.0209 

2.0 

3.0 
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Table 23. COVERAGE RATIO (MIXED EXPONENTIAL REPAIR TIME : 
ONE-SIDED, N = 100) 


Texel pes [me | moe 
* [ee [| ees [nese [ee 
snmp fam sean | sae 
* [ee [[aeon [atone 
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Table 24. LENGTH OF C.l (MIXED EXPONENTIAL REPAIR TIME 
ONE-SIDED, N = 100) 


0.099 0.0604 0.0597 0.0609 
0.2 

0.0171 0.0096 0.0098 0.0097 
0.5 

0.0168 0.0103 0.0108 0.0104 

0.1164 0.0647 0.0649 0.0642 
1.0 

0.0159 0.0136 0.0135 0.0126 
2.0 
3.0 
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APPENDIX B. TRUE AVAILABILITY TABLE 


Table 25. TRUE AVAILABILITY OF SYSTEM AT A FINITE TIME T 
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APPENDIX C. ESTIMATED AVAILABILITY TABLES 


Table 26. ESTIMATED AVAILABILITY OF SYSTEM (EXPONENTIAL RE- 
PAIR TIME: N = 50) 


SRMSE : square root of mean square error 


; SYS : 

Wass 0.8506 0.8506 0.8511 

SEM 
0.7396 0.7444 0.7455 0.7458 

0.5 
SRMSE 0.042 0.0418 0.0418 
Syl ole 0.0412 0.0437 0.0439 
0.6667 0.6704 0.6659 0.6663 

2.0 
SRMSE 0.0498 0.0417 0.0466 0.0467 
0.6665 0.6683 0.6646 0.6647 

5.0 
SRMSE 0.0469 0.0418 0.0451 0.0452 
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Table 27. ESTIMATED AVAILABILITY OF SYSTEM(EXPONENTIAL RE- 
PAIR TIME: N = 100) 


SRMSE : square root of mean square error 


SYS 
0.2 

SRMSE 0.0335 0.0272 0.0272 0.0282 
0.5 

SRMSE 0.0365 0.0313 0.0307 0.031 

0.6834 0.6831 0.684 0.6842 

1.0 

SRMSE 0.0387 0.031 0.0308 0.031 
2.0 

SRMSE ee. | 0.0318 0.0318 0.0319 
3.0 












Table 28. ESTIMATED AVAILABILITY OF SYSTEM (GAMMA REPAIR 
TIME: N = 50) 


SRMSE : square root of mean square error 


; SiS ; 
0.8346 0.8319 0.8302 0.8307 
0.2 
SRMSE 0.0485 0.0422 0.0423 0.0436 
OF tt O7232 0.7250 07238 
0.5 
SRMSE 0.0562 0.0433 0.0463 0.0454 
0.6669 0.6779 0.6775 0.6780 
SRMSE 0.0507 0.0375 0.0415 0.042 
0.6656 0.6664 0.6629 0.6630 
2.0 
SRMSE 0.0457 0.035] 0.039] 0.0392 
0.6655 0.6656 0.6634 0.6635 
3.0 
SRMSE 0.0407 0.0351 0.0379 0.0379 
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Table 29. ESTIMATED AVAILABILITY OF SYSTEM (GAMMA REPAIR 
TIME: N = 100) 


SRMSE : square root of mean square error 


SYS 
0.8337 0.8292 0.8292 0.8297 
0.2 
SRMSE 0.0363 0.0296 0.0296 0.0308 
0:5 
SRMSE 0.0387 0.0316 0.0315 0.031 
0.6672 0.6735 0.674] 0.6741 
1.0 
SRMSE 0.0388 0.0287 0.0285 0.0291 
ZO 
SRISE — 0.028 0.0282 0.0282 
3.0 
SRMSE on 0.0282 0.0276 0.0276 













Table 30. ESTIMATED AVAILABILITY OF SYSTEM (MIXED EXP REPAIR 
TIME: N = 50) 


SRMSE : square root of mean square error 


? SS 3 . 
Ue 
SRMSE 0.0319 0.0224 0.0224 0.0219 
0.8846 0.884 0.8837 0.8858 
0.5 
SRMSE 0.034 0.0222 0.0232 0.0228 
0.8738 0.8734 0.8712 0.8726 
SRMSE 0.0319 0.0245 0.0256 0.0251 
2.0 
SRMSE 0.0318 0.0261 0.0283 0.0282 
3.0 
SRMSE 0.0301 0.0265 0.029 0.029 
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Table 31. ESTIMATED AVAILABILITY OF SYSTEM (MIXED EXP REPAIR 
TIME: N = 100) 


SRMSE : square root of mean square error 


; SYS 

0.9094 0.907 0.9077 0.9081 

0.2 
SRMSE 0.0231 0.0169 0.0165 0.0164 

0.5 
SRMSE 0.023 0.018 0.0174 0.0171 
0.8752 0.8709 0.8708 0.8723 
SRMSE 0.0237 0.0196 0.0195 

2.0 

3.0 
SRMSE ~~ 2 0.0206 0.0197 0.0196 
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